Random pinning limits the size of membrane adhesion domains 
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Theoretical models describing specific adhesion of membranes predict (for certain parameters) 
a macroscopic phase separation of bonds into adhesion domains. We show that this behavior is 
fundamentally altered if the membrane is pinned randomly due to, e.g., proteins that anchor the 
membrane to the cytoskeleton. Perturbations which locally restrict membrane height fluctuations 
induce quenched disorder of the random-field type. This rigorously prevents the formation of macro- 
scopic adhesion domains following the Imry-Ma argument [Y. Imry and S. K. Ma, Phys. Rev. Lett. 
35, 1399 (1975)]. Our prediction of random-field disorder follows from analytical calculations, and 
is strikingly confirmed in large-scale Monte Carlo simulations. These simulations are based on an 
efficient composite Monte Carlo move, whereby membrane height and bond degrees of freedom are 
updated simultaneously in a single move. The application of this move should prove rewarding for 
other systems also. 

PACS numbers: 87.16.A-,87.17.Rt,75.10.Hk 



I. INTRODUCTION 

The fate of living cells is regulated through interac- 
tions with other cells and with the extracellular matrix 
(ECM). Through receptor-Hgand bonds formed by spe- 
cific proteins the cell adheres to the ECM forming ad- 
hesion domains (clusters of closed bonds). Not all ad- 
hesion domains are focal adhesions but these are par- 
ticularly well studied and relevant. Focal adhesions are 
involved in the transmission of signals and mechanical 
forces, and play key roles in cell anchorage and migra- 
tion [2 [2]. Consequently, understanding how adhesion 
domains form, and the factors that control their size, 
shape, and growth [3l |4] , is of profound practical impor- 
tance. Extensive studies have been performed on the- 
oretical models for single bond dynamics [S], collective 
dynamics of discrete bonds [6Hl0]. and employing effec- 
tive potentials pdlll3j : as well as experimentally on cell- 
mimetic model systems [TH [TS] such as lipid bilayer vesi- 
cles with embedded ligands brought in the vicinity of 
receptors tethered to supported membranes [I6l - f20] . 

While theoretical models predict a macroscopic phase 
separation, adhesion domains in cells are typically finite 
in size. It has been proposed that trapping in membrane 
corrals (or compartments) reduces the mobility of recep- 
tors [H] , which becomes a limiting factor for the growth 
of adhesion domains. The purpose of this paper is to pro- 
vide an alternative explanation based on the observation 
that, due to the interactions with the ECM, membrane 
height fluctuations are locally suppressed. These "pin- 
ning sites" induce quenched disorder of the random-field 
type which prevents macroscopic domain formation in 
d — 2 dimensions P^25j . Hence, based on the funda- 
mentals of statistical physics alone, adhesion domains of 
finite size are implied in a wide variety of contexts. The 
mere presence of random-field disorder provides a robust 
mechanism controlling the size of adhesion domains, irre- 



spective of the details of the many and complex molecular 
processes between cell and ECM. 

The generic situation that we envision is sketched 
in Fig. fl[a), which shows a fluctuating membrane adher- 
ing to a substrate via receptor-ligand bonds. The crucial 
ingredients are the pinning sites at which the membrane 
height is assumed to be fixed. In a biological cell, such 
pinning sites correspond, e.g., to the locations where the 
cytoskeleton anchors to the membrane. For an experi- 
mental verification of our predictions in vitro we propose 
a bilayer adhered to a substrate or supported membrane 
with a (small) fraction of the receptor-ligand bonds per- 
manently closed and of high stiffness [Fig. llTb)]. We will 
show that in situations resembling those of Fig. [l] adhe- 
sion domains remain finite in size for any finite concen- 
tration of pinning sites, provided the spatial distribution 
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FIG. 1: (Color online) Sketch of two typical situations in 
which the formation of macroscopic adhesion domains is pre- 
vented due to the presence of random fields: (a) A mem- 
brane of a biological cell adheres to the ECM via receptor- 
ligand bonds, which may dynamically open and close. The 
membrane is freely fiuctuating except at the anchoring sites 
(pinning sites) attaching the membrane to the cytoskeleton. 
(b) An analogous in vitro example in which the pinning sites 
are formed by stiff and permanently closed bonds (tethers) 
between membrane and substrate. 



of pinning sites is random. 



II. MODEL 

To show how the random-field disorder in the specific 
adhesion of membranes comes about we consider a class 
of simple models as reviewed in Ref. HSl These coarse- 
grained models contain the minimal ingredients that we 
require to address the influence of membrane pinning on 
the statistics of adhesion domain formation. In partic- 
ular, we assume the membrane to be a two-dimensional 
sheet characterized by a bending rigidity, and to be in 
thermal equilibrium with its surroundings. The ligand- 
receptor bonds and pinning sites are treated as point par- 
ticles. We neglect all active processes and, in the case of 
focal adhesions, the influence that stresses applied to the 
substrate may have on the development of adhesion do- 
mains "7. 

In the Monge representation a membrane patch with 
projected area A = L"^ is described through its separa- 
tion profile h(r) describing the height of the membrane 
at position r = {x,y), ^ x,y < L, measured with re- 
spect to some (arbitrarily chosen) reference height. The 
effective Hamiltonian reads 



na[h{r)] = 



d^i 



{l[^'h{r)r + l[Hrr}. (I) 



The first term governs the bending energy, with n the 
membrane bending rigidity, which we assume is the dom- 
inant contribution to the Helfrich energy ^7\. The 
second term is the lowest-order expansion of the non- 
specific interactions between membrane and substrate 
whose strength is denoted 7. In our treatment, the min- 
imum of the non-specific potential is thus taken to be 
the reference height from which h{r) is measured. These 
non-specific interactions arise due to volume exclusion, 
van der Waals forces, and the possible formation of an 
electrostatic double layer, as well as an effective pressure 
due to the restricted volume the membrane can move in. 
The parameters k and 7 define a length, ^y = {n/^y^^, 
which sets the scale over which the membrane height fiuc- 
tuations are correlated |8], 

m(|r - r'l) ^ {h{r)h{r'))o « ei"(|r - rV^n), (2) 

see appendix P] The scaling function is defined in terms 
of the Kelvin function, u{x) — — (4/7r) keio(a::); it obeys 
m(0) = 1 and decays to zero exponentially fast. The am- 
plitude of the correlations is set by the thermal roughness 
f_L = ^||/>/8/Jk, /? = {kBT)~^, with temperature T and 
Boltzmann constant fee- The brackets (• • •)o denote the 
thermal average with respect to the Hamiltonian Eq. (IT]) , 
i.e., in the absence of ligand-receptor bonds and pinning 
sites. 

Now imagine a set of TV receptor-ligand pairs at po- 
sitions R = {ri, r2, . . . , r^v} embedded in the membrane 



and substrate that can form bonds. The arguably sim- 
plest model is to assign a linear energy to closed bonds 
with stiffness a [5] . The total Hamiltonian then reads 



JV 



Hat = Ho + ^ bi [ah{ri) - Sb] , 



(3) 



where 6^ = if the ith bond is open, and 6^ = 1 if the 



bond is closed. The parameter £b 



ad is the 



shifted binding energy the system gains through form- 
ing a bond, where e^ is the bare binding energy and d 
is the separation between substrate and the minimum of 
the non-specific potential. In addition to the receptor- 
ligand pairs, we assume the presence of n pinning sites 
at positions P = {rjv+i, . • • ,J^m} where the membrane 
height is fixed to h — hp. We take the two sets R and P 
to be disjoint with their union thus holding AI — N + n 
distinct sites. 



III. THEORETICAL MAPPING 

We now calculate the free energy as a function of the 
bond variables {bi} through integrating out the height 
fiuctuations under the constraints h{ri) = hi for sites 
Ti e R and h{ri) — hp for r^ G P. We follow the 
standard procedure and implement these M constraints 
through (5- functions ^28j; for the detailed calculation see 
appendix |BJ 

For clarity of the presentation, in the following we con- 
sider the case hp = 0. The free energy then reads 

G({6J; RUP) = -/3-1 In / d/ii • • • dh^ x 



exp 






)ij f^if^j 



N 

I3^bi[ahi -Eb] 

i=l 



(4) 



The dependence on the positions R and P is encoded 
in the matrix m, whose components are given by rriij = 
m{\ri — Vj\). Performing the final integrations we obtain 



G^ 



N N 

EE 

1=1 j=i+l 



N 



Eijbibj -'^Hibi, Eij = Pa^rfiij, (5) 



with ^.i = eh+l3a'^mii/2, and where we have introduced a 
new matrix rh. While Eq. Q contains the full matrix m, 
the sums run only over the first N sites corresponding to 
the bonds and excluding the pinning sites. The matrix m 
is obtained by inverting the submatrix formed by the first 
N rows and N columns of m~^. Note also that we have 
dropped an additional term in Eq. (15 ) which depends on 
the geometry of both the bonds and the pinning sites but 
not on the state {bi} of the bonds. 

The key result here is that the free energy in Eq. (Isl) is 
isomorphic to the Ising lattice gas with couplings Sij ex a^ 
and a site- dependent effective chemical potential ^i that 



also depends on the stiffness a. There are two effects due 
to membrane undulations: First, single bond formation is 
assisted (^^ > eb) since the system can access configura- 
tions with lower energy. Second, bonds couple in a man- 
ner that enhances clustering: a closed bond pulls down 
the membrane locally making it easier for nearby bonds 
to also close. In such bound patches (adhesion domains) 
membrane fluctuations are hindered and, therefore, the 
entropy is decreased. The phase behavior of the system 
is thus determined by the competition between this loss 
of entropy, the mixing entropy of bonds, and the gain in 
binding energy. For small ^ii only nearest neighbors inter- 
act directly; by increasing ^y (i.e., for stiffer membranes) 
the thermal roughness ^_l determining the strength of the 
coupling is diminished but more and more bonds become 
coupled. 
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FIG. 2: (Color online) Graphical aid to identify the non- 
vanishing components of the matrices A^ and A^ appearing 
in Eq. (pi with bond sites i and j: (a) For i = j the sites 
fci , . . . , fc4 that are nearest neighbors of i are counted, (b) For 
the case that i and j are next-nearest neighbors the sites fcs 
and fcg that are common nearest neighbors of both i and j 
are counted. 



A. Membrane without pinning 

The result Eq. ^ holds for any geometry of bonds 
and pinning sites. In the absence of pinning sites (n — 
0) clearly rfi = m and the chemical potential becomes 
spatially uniform, /x^ = £b + eo/2 = fJ-o, where 



£o 



/3(«a)' 



aV(8v^) 



(6) 



is the effective coupling energy due to the membrane un- 
dulations. We now specialize to the situation where the 
positions R of receptor-ligand pairs form a regular square 
lattice with lattice spacing a ^ ^|| . The phase behavior 
in this case is known [6l|TT]: For sufficiently high a > a* 
there is a first order phase transition from a bound state 
{(j> ^ I) with a high density of closed bonds 



1 " 



(7) 



to an unbound state with a low density of closed bonds 
{(f) ~ 0). Precisely dX a = a* the system shows criti- 
cal behavior which, by virtue of the mapping of Eq. ([5| , 
belongs to the universality class of the d = 2 Ising model. 



B. The pinned membrane 

We now come to the main result of this paper, where 
the fate of the mapping of Eq. ([s]) in the presence of pin- 
ning sites {n > 0) is considered. In this case, we must dis- 
tinguish between the set R of sites where receptor-ligand 
bonds may form, and the set P of sites at which the mem- 
brane is pinned. We furthermore restrict our calculation 
of m to nearest and next-nearest neighbor interactions. 
We first write m = ^]^ [1 -h ca A + €bB], where the M x M 
matrices A and B correspond to nearest and next-nearest 
neighbor interactions, respectively. The components Aij 
(Bij) equal one if the two sites i and j are nearest (next- 
nearest) neighbors, and zero otherwise. The coefficients 



are set by the value of the scaling function at the nearest 
and next-nearest neighbor distance, ea = u{a/M and 
eb = ■'^( v2a/'?||), respectively. We expand the inverse as 



(m/ei)- 



l-EAA-eBB + eiA^ 



(8) 



where all terms with coefficients < e\ have been dropped. 
We now replace, in Eq. (|8|, the M x M matrices {1, A, 8} 
by the N x N sub-matrices {1,A, B} obtained by dis- 
carding the n upper rows and columns corresponding to 
pinning sites; taking the inverse of the resulting matrix 
yields the desired matrix rti that appears in the mapping 
of Eq. ^ 



rh/ei « i + eaA + ebB + el ( A^ - A^ 



(9) 



where all higher order terms have again been dropped. 

Due to the pinning sites, the last term in Eq. (|9| does 
not vanish since the "hat" and "square" operations do 

not commute. First we evaluate {l^'^)ij — '^k=i ^ik^kj 
where, in the summation, only two terms survive (Fig.^. 
The first term corresponds to i = j; a non-zero contri- 
bution implies that k needs to be a nearest neighbor of 
this site. On a square lattice there are four such sites 
labeled fci , . . . , ^4 in Fig. l2](a) . The second term arises 
when i and j are distinct. The only combination with a 
non-zero contribution is when i and j are both nearest- 
neighbor of the same site fc, which implies that i and j are 
next-nearest neighbors. Two such sites can be identified, 
labeled k^ and fcg in Fig. [2jb). Hence, 



(A^). 



4(5, 



2B, 



(10) 



with 5ij the Kronecker symbol. The components of 
the square of the reduced nearest-neighbor matrix are 
(A^)ij = X]fc=i ^jfc^fcj J where the sum over k now ex- 
cludes the pinning sites. We can identify the non- 
vanishing terms as in Fig. [2] provided we ignore the sites 
ki that are pinned. We thus obtain 



(A^).. 



i'i-ndSij + XijBi, 



(11) 



with stochastic variables rji and Xij set by the local envi- 
ronment of pinning sites. In particular, rji is the number 
of nearest-neighbors of site i that are pinned; the possi- 
ble values of Xij = 0, 1, 2 correspond to, respectively, the 
case where both k^ and fcg are pinned, only one of those 
sites is pinned, and neither one of them being pinned. 

We now have all the ingredients needed to discuss the 
mapping of Eq. (l5| in the presence of pinning sites. The 
first observation is that the nearest-neighbor coupling is 
not affected. For sites i and j that are nearest-neighbors, 
Eij = SoeA as before, by virtue of Eq. Q. In contrast, 
the chemical potential is affected, and now depends on 
the local environment via r]i 



IJ-i = Mo - ??iA, A = £061/2, 



(12) 



where mq is the effective chemical potential in the absence 
of pinning sites. Provided the pinning sites are immo- 
bile and randomly distributed, Eq. (12) corresponds to 
a quenched random-field. The effective chemical poten- 
tial is thus reduced, implying that the closing of bonds 
has become more difficult. The physical picture is that, 
since the membrane is pinned to the minimum of the non- 
specific potential (/ip = 0), closed bonds cannot "pull 
down" the membrane as easily as before. This impedes 
the above mentioned facilitation of bond formation due 
to undulations, and consequently the effective chemical 
potential is reduced. Choosing a sufficiently negative 
/ip < 0, the opposite situation of an increased local chem- 
ical potential may also be realized. 

We typically consider a low pinning density p = 
n/N ^ 1, such that rji effectively becomes a binary ran- 
dom variable with values rji — 0, 1. In this limit, the aver- 
age chemical potential is [/i^] = /io^ pA, where [•] denotes 
a disorder average (i.e., an average over many different 
samples of pinning sites). The random- field strength is 

set by the disorder fluctuations \/[m?] ~ [m*]^ ~ y/~P^i 
which thus is weak compared to the nearest-neighbor cou- 
pling. Nevertheless, given that an infinitesimally weak 
random-field is sufficient to destroy macroscopic domain 
formation in two dimensions |22H25) , we expect that even 
a low pinning density will drastically affect adhesion do- 
main formation. Finally, note that in addition to the 
dominant random-field disorder, the pinning sites also 
induce a marginal perturbation. For next-nearest neigh- 
bors i and j, the coupling e^ = £o[eB + ^\{Xij ~ 2)] also 
becomes a random variable corresponding to random- 
bond disorder, which does not destroy macroscopic do- 
main formation. Possible interactions between receptor- 
ligand bonds, e.g., due to size mismatch, will change the 
couplings Eij but not the fact that a random-field is in- 
duced. 



IV. SIMULATIONS 
A. Model and methods 

1. Discretized membrane model with pinning sites 

We now perform Monte Carlo (MC) simulations us- 
ing a discretized version of our model Hamiltonian. We 
mostly simulate on periodic N = L x L lattices, and to 
each lattice site i we assign a real number hi to denote 
the local membrane height, and a bond variable bi to de- 
note whether the bond at the site is open {bi — 0) or 
closed (bi = 1). The Hamiltonian may then be written 
as 



n 



N 



s:[i<-= 



hif + ^h^, + b,{ah, - eh)\ , (13) 



where the sum extents over all lattice sites. To compute 
the Laplacian, we use the finite-difference expression 



V^/i., = Xi = -Ahi 



jenn(j) 



(14) 



where the sum in the last term is over the (z — 4) nearest 
neighboring (nn) sites of site i. In this section the inverse 
temperature /3 is absorbed into the coupling constants 
{k,7, Q:,eb} of Eq. (13), while the lattice constant will 



be the unit of length. 

To incorporate pinning, a fraction p of randomly se- 
lected sites is placed in the set P of pinning sites. The 
sites i G P have their height variable set to hi = hp 
at the start of the simulation, and these heights are not 
permitted to change during the course of the simulation 
(i.e., MC moves that change the membrane height are 
not applied to the pinning sites). In contrast to the the- 
oretical derivation we do allow for bonds to open and 
close at the pinning sites. In the simulations, the sets R 
and P therefore overlap. This does not affect the phase 
behavior of Eq. ( [l3| but makes the data analysis easier 
since the available area for bonds then always equals L^ 
as opposed to {1 — p)L^. 



2. Composite Monte Carlo move 

The "standard approach" to simulate Eq. ( jTs] ) is to 
use MC moves that either (i) propose a small random 
change to a randomly selected height variable, or (ii) 
change the state of a randomly selected bond variable, 
and to accept these changes with the Metropolis crite- 
rion [6^ . This approach is not efficient because the height 
and bond degrees of freedom are correlated: at sites con- 
taining a closed bond, the membrane height will be lower, 
and vice versa. Hence, after a proposed change in bi, the 
corresponding height hi is likely to be energetically un- 
favorable, and so move (ii) has a high chance of being 
rejected. 



To circumvent this problem, we use a composite MC 
move whereby the bond and height degrees of freedom 
are changed simultaneously. The key observation is that 
Eq. ( 13 1 is quadratic in the height variables. For a given 
site i, imagine to replace the corresponding membrane 
height by hi i^ h[ ~ hi + 6. The energy as function 
of 5 is quadratic, E{S) — A2S^ + AiS + constant, with 
coefficients given by 



A2 = 10k + 7/2, 



Ai — ^hi + abi — K 4Ai 



jenn(i) 



(15) 



(16) 



Hence, there is an optimal deviation, Jopt = — A1/2A2, at 
which the energy becomes minimized. In our MC simula- 
tions, we exploit this property by selecting the membrane 
height deviations 5 from a Gaussian distribution around 
the optimal value 



P{5\5o 



Vi) 



-(5-<5opt)V2<T" 



(17) 



where standard deviation a^ — I/2A2 is used (other 
choices for a^ are valid also, but we believe this one is 
optimal as it closely matches the thermal fluctuations). 

The composite MC move that we use to change the 
state of bond variables proceeds as follows: 

1. Randomly select a lattice site i, and compute the 
optimal height deviation (5opt. 

2. Change the state of the bond variable bi 1— > 6J, 
and compute the new optimal height deviation 5'^. . 
Propose a new height, hi ^-^ h^ = K + S, with S 
drawn from P{5\6'^p^) of Eq. ([l7|. 



3. Accept the proposed 
Metropolis criterion 



values {b'ijh'i) with the 



Parr. = min 



. P{0\Sopt) . 

' P{SK,t) 



(U'^-Un) 



(18) 



with 'Hat the energy of the configuration at the start 
of the MC move, and T-L'j^ the energy of the pro- 
posed configuration (the ratio of Gaussian proba- 
bilities is needed to restore detailed balance). 

Note that our composite move is ergodic, and thus by 
itself constitutes a valid MC algorithm. Nevertheless, we 
still found it useful to also implement the non-composite 
variant, whereby the membrane height is updated with- 
out changing the corresponding bond variable. In terms 
of the composite move above, this corresponds to 5'^ = 
<5opt , while in step (2) the bond variable bi is not changed. 
In what follows, composite to non-composite moves are 
attempted in a ratio 1:2, respectively. 

In the MC moves above, we restrict the selection of the 
optimal height value hi to the single site i. An obvious 
generalization is to also optimally select the height val- 
ues of nearby sites, i.e., on a Z x Z plaquette around site i. 



The above moves correspond to Z = 1, but we have used 
the version with / = 3 also; the latter is slightly more 
efficient in cases where bound and unbound membrane 
patches coexist. Obviously, the case Z > 1 is more com- 
plex to implement, as it involves minimizing a quadratic 
form of P variables. Furthermore, the accept criterion 



Eq. ( 18 1 needs to be modified (the prefactor now becomes 
the product of P Gaussian probability ratios). 



3. Order parameter distribution 

A key output of our simulations is the order param- 
eter distribution P{(t>) (OPD), which is defined as the 
probability to observe a system with a fraction of closed 
bonds 0, with (j) given by Eq. ([7]). We emphasize that 
P{<t>) depends on all the coupling constants that appear 
in Eq. (13), as well as on the system size L. To ensure 
that P(^) is sampled over the entire range < ^ < 1, 
we combine our simulations with an umbrella sampling 
scheme [33]. We also use histogram reweighting ^30] in 
the binding energy: having measured P{(j)) for some value 



Eb 



£ 



we extrapolate to different values Eb ~ £b using 



the relation 



•\£h 



zl) (X P{ 



'\eh 



eO)ei^(4-eS)0. (19) 



In a similar way we also use histogram reweighting in the 
coupling constant a, which is slightly more complex to 
implement as it requires separate storage of the fluctua- 
tions in ^^ bihi, i.e., the third term of Eq. ( 13 ). 

We emphasize that, in the presence of pinning, P{(j)) 
may also depend on the particular sample of pinning 
sites. For an accurate analysis, it then becomes neces- 
sary to average simulation results over many different 
random positions (samples) of the pinning sites. 



B. Membrane without pinning 

We first simulate a membrane without pinning {p = 0) 
using K = 7 = 1 in the Hamiltonian of Eq. (Il3|. This 
case was considered extensively in Ref. W^ the main con- 
clusion being that macroscopic adhesion domains are ob- 
served for a> a* . We revisit this case to also determine 
the universality class, as well as the line tension between 
coexisting domains. 



1. Critical behavior 

We first determine the critical value a* via finite-size 
scaling of the order parameter m = (|(5(/'|), the suscepti- 
bility X = L^ {{^4''^} — (l'50l)^)i and the Binder cumulant 
U4 = {S(t)'^f/{S(j)y. Here, Scj) = (j) - {<f>), with thermal 
averages (</)") = /^ 4''^P{<j))d(f), where it is assumed that 
P{4i) is normalized. We emphasize that m, x, U^ are to be 
computed at the coexistence value of the binding energy 
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2. Symmetry line 

In Fig. IsFd), the variation of eP with a as obtained 
in our simulations using Eq. ( |20[ ) is shown. Coexis- 
tence of the bound and the unbound phases occurs along 
this symmetry line, Eh = £'^{a), which implies that the 
Hamiltonian is invariant under "swapping" the two coex- 



isting phases. For the discrete Hamiltonian Eq. ( 13 ) this 
corresponds to the operation 



{S„bi) ^ {-S,,l-bi), 



(21) 



where 6i = hi ~ (h) is the height deviation at lattice 
site i around the mean membrane height (h) . A straight- 
forward calculation shows that, in order for "H^r to be 
invariant, we are left with the condition 



1 



N 



(2(/i)+a)-^<5, = (aW-en(l-20). 



(22) 



FIG. 3: (Color online) Finite-size scaling results for the mem- 
brane without pinning sites. The scaling plots of the order 
parameter (a) and susceptibility (b) confirm 2D Ising univer- 
sality, as indicated by the collapse of the data from different 
system sizes. The intersection of the Binder cumulant curves 
for different system sizes (c) yields an estimate of a* consis- 
tent with the scaling plots. The variation of e^ with a (d) 
conforms to the theoretical prediction, where the dot marks 
the critical point. 



^b 



An accurate numerical criterion to determine 



the latter is to tune eb such that the fluctuations in 
become maximized 1311 



') - {<Pf 



max, 



(20) 



which may conveniently be done using the histogram 
reweighting formula of Eq. (19 1. 



In Fig. Isla), we show the scaling plot of the or- 
der parameter [32], i.e., curves of mL^^^ versus tV-/'^ ^ 
t — a/a* — 1, using 2D Ising critical exponents (3 — 1/8, 
V = 1^ and with a* — 2.665 obtained by tuning until the 
curves for different L collapsed (note: we use the "stan- 
dard symbol" j3 for the order parameter critical exponent, 
which is not to be confused with the inverse tempera- 
ture). The fact that the data collapse confirms 2D Ising 
universality, and our estimate of a* is in good agreement 
with Ref. W. In Fig. Isjb), we show the corresponding 
scaling plot of the susceptibility, using the 2D Ising value 
7 = 7/4, while for a* the above estimate was used. A 
data collapse is again observed, providing further con- 
firmation of 2D Ising universality. In Fig. IsFc) , we plot 
the Binder cumulant [/4 versus a. In agreement with a 
critical point, the curves for different L intersect at a* . 



Hence, the mean height along the symmetry line obeys 
{h) — —a/2. Comparing this result with an alterna- 
tive calculation in appendix O we find that (</>) = 1/2 
along the symmetry line, as expected. Moreover, from 
the second condition we obtain the symmetry line: e'^ — 
—a^/2. As Fig. p[d) shows, the agreement with the sim- 
ulation result is excellent. 

In the effective model Eq. (|5|, spin-reversal symmetry 
corresponds to the operation 6^ i— > 1 — 6^. The symmetry 
line is now determined through 



^^'i^ = 



\^' 



= £h 



£o/2, 



(23) 



JT^i 



where the sum runs over all sites j excluding i. Plug- 
ging in the definition Eq. ([6]) for the coupling energy Eq, 
the binding energy at coexistence is e^ = —ca^ with 
c = (^^/2)(1 -I- 4eA + 4eB + ■■■)• The prefactor c appar- 
ently depends on the correlation length ^n determining 
the interaction range, but should nevertheless converge 
to c = 1/2. We have checked for i^|| = 1 that this is 
indeed the case [Fig. Isld)]. 

At the coexistence binding energy, e'^ = — a^/2, the 
order parameter probability P(0) is symmetric about 
(j) = 1/2 [Fig. Ilia)]. The latter reflects the spin-reversal 
symmetry of the Ising model, which persists in the mem- 
brane model as becomes evident from the theoretical 
mapping. We emphasize that this applies to the mem- 
brane without pinning sites: in the presence of pinning 
sites, spin-reversal symmetry is generally broken. Instead 
of using the density of closed bonds (j) as the order pa- 
rameter, we could also have used the membrane height 
per site h = [l/V^) ^^ hi [6| since the latter is directly 
coupled to the density of closed bonds. This can also be 
seen in Fig. llFc), where we show the same snapshot as in 
(b) , but this time color-coded according to the membrane 
height. Furthermore, in a canonical (fixed (j)) simulation, 
and using the symmetry value (p = 1/2, the binding en- 
ergy always assumes the coexistence value e'^ (for the 
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FIG. 4: (Color online) (a) Free energy In P{<j)) at coexis- 
tence for two values a S> a* using a rectangular 50 x 100 
simulation lattice. The vertical double-arrow indicates the 
free energy cost AF of interface formation. Note the pro- 
nounced flat region between the peaks, (b) Snapshot of the 
membrane at = 1/2 and a = 2.9, where the dots indicate 
closed bonds, (c) The same snapshot, but color-coded ac- 
cording to the membrane height. The snapshots (b,c) clearly 
indicate that, in regions with a high density of closed bonds, 
the membrane is "pulled down". 



Applying this equation to the distributions of Fig. |4][a) , 
we obtain a ~ 0.34(0.48) for a = 2.9(3.0), in units of 
kgT per lattice spacing. As expected, a decreases upon 
lowering a. 



4- Summary 

For the membrane without pinning, our simulation re- 
sults are thus fully consistent with the theoretical pre- 
diction that such a system should map onto the 2D Ising 
lattice gas. By means of finite-size scaling, 2D Ising crit- 
ical exponents are confirmed, the quadratic variation of 
the coexistence binding energy with a is recovered, and 
we observe the analogue of spin-reversal symmetry. Fur- 
thermore, our estimate of a* is in good agreement with 
Ref. !i6 , providing an important consistency check for the 
composite MC moves proposed in this work. Finally, for 
a 3> a*, we observe a first-order phase transition be- 
tween bound and unbound phases, with an associated 
coexistence region where macroscopic domain formation 
occurs. 



Ising model with conserved order parameter |32j , the ana- 
logue of this condition is that, at zero magnetization, the 
external field is zero). In a grand-canonical simulation, 
i.e., where (j) is allowed to fluctuate, and using the coexis- 
tence binding energy e'^, the probability distribution of 
the membrane height P{h) will thus be symmetric about 
h = -a/2. 



3. Line tension 

Next, we consider a 3> a*, where the transition is 
strongly first-order. In Fig.l4|a), we show lnP(0) at co- 
existence for two values of a, using a rectangular 50 x 100 
simulation box (note that lnP((/)) may be regarded as 
minus the free energy of the system). We observe that 
lnP(0) is profoundly bimodal, which indicates two-phase 
coexistence. The left peak corresponds to the unbound 
phase (low density of closed bonds), the right peak to 
the bound phase (high density of closed bonds), while 
in the region "between the peaks" both phases appear 
simultaneously. The latter follows strikingly from snap- 
shots, see Fig. |4]^b), which was taken at ^ = 0.5. For 
this value of <f), the phases arrange in two slabs, with two 
interfaces perpendicular to the longest edge, as this min- 
imizes the total length of the line interface (whose length 
then equals Ltot = 2_Lshort, with Lghort = 50 the short- 
est edge of the simulation box) . Provided the interfaces 
do not interact, there is a region over which (p can be 
varied without any free energy cost. This, apparently, is 
the case here, as the distributions lnP(0) are essentially 
flat in their center regions. Following Binder |33j . we 
may then relate the free energy barrier AF, indicated in 
Fig. l4[a) for a = 3.0, to the line tension a = AF/Ltot- 



C. The pinned membrane 

We now consider a membrane with a finite concentra- 
tion of pinning sites p > 0. We continue to use k = 7 = 1 
in Eq. ( |13[ ), for which the case without pinning was just 
described (Section IVB ). In what follows, wc restrict our- 
selves to p ^ 1, i.e., the limit of low pinning concentra- 
tion. We believe this to be the biologically most relevant 
situation, as well as the physically most interesting one 
(in the opposite limit p -^ 1, the pinning sites would com- 
pletely freeze the membrane, thereby trivially preventing 
any membrane-mediated phenomena from taking place) . 



1. Numerical evidence for random-field disorder 

The key prediction of the theoretical mapping, Eq. ([5| , 
is that, in the presence of pinning sites, the chemical po- 
tential becomes a quenched random variable (i.e., depen- 
dent on the spatial location in the sample, but without 
any time dependence). We thus expect "special regions" 
in the membrane where bonds prefer to close, which are 
those regions where the local chemical potential excess is 
positive. To test whether such regions can be found, we 
perform canonical MC simulations at a fixed fraction of 
closed bonds (j) — 1/2. For each lattice site i, we measure 
the thermally averaged bond occupation variable (&i), 
and use these to compute the spatial fluctuation [3T| l34] 
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where the sum is over all lattice sites. In the absence of 
spatial preference, {hi) = 1/2 for all sites, implying that 
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FIG. 5: (Color online) (a) The "running average" of xs, 
Eq. 1 24 1, which measures the spatial fluctuation in the ther- 
mally averaged bond occupation variables (bi). The solid 
curves (blue) correspond to a pinned membrane for a = 
2.0, 2.5, a*, 3.0, 3.5,4.0 (from bottom to top) where a* is the 
critical point of the unpinned membrane. As a increases, 
Xs — >■ 1/2 (horizontal line). The dashed curve (red) shows 
Xs for the membrane without pinning at a = 2.5, in which 
case Xs decays to zero. (b,c) "Color-maps" of the thermally 
averaged bond occupation variables (bi) for a = 2.0 (b), and 
a = 4.0 (c). In both cases, the same pinning configuration 
(red dots) is used. (All data in this figure refer to L = 100; 
results for the pinned membrane use p = 0.05 and hp — 0). 



Xs — 0. However, if sites in certain regions prefer closed 
bonds, then (bi) will be distinctly different from 1/2, and 
consequently xs > 0. 

In Fig. pi a), we show "running averages" of xs a-s a 
function of the number of MC moves r. For the un- 
pinned membrane, xs decays to zero, since here there is 
no spatial preference. In the absence of spatial prefer- 
ence, Xs reflects the statistical error of our simulation, 
and therefore decays oc t~^/^. In contrast, for the mem- 
brane with pinning sites, xs saturates to a finite plateau 
value. Note that this happens for all values of a con- 
sidered, including those values below the critical point 
a* of the unpinned membrane. By increasing a, the 
plateau value increases, which means that the spatial 
preference becomes stronger. When a is small, thermal 
fluctuations still permit "excursions" of bonds into re- 
gions where the local chemical potential is unfavorable. 
In this case, (bi) does not deviate much from 1/2. As a 
increases, these thermal fluctuations are "frozen out" . In 
the limit a — > oo, the bond occupation variables are set 
by the groundstate of Eq. (13| (that is: for a given config- 
uration of pinning sites, the variables bi are determined 
by energy minimization). In this limit, (bi) is either 
or 1, and the fluctuation xs — 1/2- As can be seen in 
Fig. p[a) , this limiting value is indeed approached as a 
increases. 

However, already for a much smaller, the groundstate 
reveals itself. To demonstrate this, we show "color-maps" 




FIG. 6: (Color online) The (normalized) correlation function 
g(r) of the local chemical potential excess A/ii, for several 
values of a. The correlations quickly decay to zero, which 
confirms the prediction of the mapping that A/ii is a spatially 
random variable (results refer to the pinned membrane with 
p = 0.05, hp = 0, and L = 100). 



of the thermally averaged values {bt) |3ll |34j, for two 
values of a, but using the same conflguration of pinning 
sites with /ip = [Fig. [5Fb,c)]. Clearly visible is that 
the adhesion domains (dark regions) prefer the same lo- 
cations for both values of a. The only difference is that, 
with increasing a, the preference becomes stronger. The 
structure of Fig. [sjjc) is already close to the groundstate, 
since the majority of {bi) are already close to either or 
1. But also for the smaller value of a, the groundstate 
is visible, albeit somewhat "blurred". Note that a in 
Fig. p[b) is considerably below a* of the unpinned mem- 
brane, i.e., the presence of the groundstate persists deep 
into the high-temperature region. 

Our finding that, in the presence of pinning sites, {bi) 
generally deviates from 1/2 may be conceived as a local 
chemical potential excess (field) at site i given by 



A^i = In 



{h 



i-ih)' 



(25) 



where the canonical ensemble with cj) = 1/2 is assumed. 
Furthermore, the "color-maps" of Fig.[5Fb,c) suggest that 
A/ii is a spatially random variable. Since the chemical 
potential is isomorphic to an external field in the Ising 
model, this indeed corresponds to random-field disorder. 
Consequently, the pinned membrane belongs to the uni- 
versality class of the two-dimensional random-field Ising 
model. This implies that the "freezing out" of the ther- 
mal fluctuations with increasing a in Fig. [5] is a gradual 
process, i.e., there is no phase transition associated with 
it (and nothing special happens at the critical point a* 
of unpinned membrane). As is well known, the random- 
fleld Ising model in d = 2 dimensions does not support 
any phase transition [35j. 



To confirm that A/z^ truly is a spatially random vari- 
able, we have measured the (circularly averaged) correla- 
tion function g{r) ex [A/ZjA/Uj]', where [•]' denotes a spa- 
tial average over pairs of sites {i,j) a distance r apart. 
As can be seen in Fig. pi g{r) quickly decays to zero. For 
randomly distributed pinning sites, the spatial correla- 
tions in the local chemical potential excess are thus short- 
ranged. Note also that g{r) depends only weakly on a. 
This is consistent with the "color-maps" of Fig. [5jb,c), 
whose overall topology (i.e., the shape of the regions) is 
also remarkably insensitive to a. The insensitivity to a 
shows that the sign of the local chemical potential excess 
in Eq. (25 1 is determined exclusively by the properties of 



the membrane and the pinning sites (most notably the 
positions of the latter, and the pinning height hp). The 
picture that one should have in mind, therefore, is that 
of ligand-receptor bonds "diffusing" through a chemical 
potential landscape, but the landscape itself is quenched, 
i.e., the bonds do not modify nor shape it. Incidentally, 
for hp = 0, Fig. [5]^b,c) shows that the local chemical po- 
tential excess is such that the pinning sites "repel" closed 
bonds. By using a sufficiently negative pinning height, 
the reverse situation can be realized also. 



2. Adhesion domains are finite: Imry-Ma argument 

For the pinned membrane, Figs. [5] and [6] clearly show 
that the local chemical potential excess is a quenched 
random variable, thereby confirming the theoretical pre- 
diction. This rigorously rules out the formation of large 
(macroscopic) adhesion domains, by virtue of the Imry- 
Ma argument P^25| . To see this, consider a cluster 
of closed bonds of linear size /, thereby containing ex l'^ 
bonds, and imagine to insert the cluster at some location 
in the sample. The average chemical potential excess 
(per site) over the cluster is zero, but with fluctuations 
that decay conform the central limit theorem 



1 



J2 AA., = o±c//'^/^ 

i^cluster 



(26) 



with A/ii given by Eq. (25), and C a constant. Hence, 



there exist "preferred regions" where the cost of insertion 
is reduced by an amount ex ^'^'^. In d = 2 dimensions 
(but not in d = 3 |35 1 136 ) 1. this is sufficient to compen- 
sate the line tension, which scales ^ l'^~^ . In the presence 
of pinning sites, adhesion domains thus no longer strive 
to minimize the length of their line interface. Instead, 
they seek out those regions in the sample where the lo- 
cal chemical potential excess is most favorable, precisely 
what is observed in Fig. [5jb,c). We thus obtain a stable 
multi-domain structure. Note the sharp contrast to the 
case without pinning, where adhesion domains do min- 
imize their interface, thereby growing macroscopically 
large [Fig. [i] . 



3. Domain size statistics 

Having argued that adhesion domains in the presence 
of pinning sites remain finite, we now present a more 
quantitative analysis of their size. As the domain struc- 
ture is essentially fixed by the pinning sites, and much less 
by thermal effects, we restrict ourselves to a = 3.5, and 
vary only (i) the pinning concentration p, and (ii) the lat- 
eral correlation length ^y of the membrane fluctuations. 
The simulations in this section are again performed in 
the canonical ensemble {(j) = 1/2), and pinning height 
hp = {) (pinning sites thus "repel" closed bonds). The 
remaining fixed parameters are 7 = 1 for the non-specific 
potential, and system size L = 100. 

For a given sample i — 1, . . . , X of pinning sites, we 
generate a series t = 1, . . . , T of equilibrated snapshots. 
For each snapshot, we compute the typical domain size 
R^^^ = 2nJSt4k)dk/JkStAik)dk [32], where St^k) is 
the (circularly averaged) static structure factor; disorder 
[•] and thermal (•) averages are then computed as 



[{Rl 



1 ^ /i ^ V 

i=l \ t=l / 



(27) 



We will primarily be concerned with the average domain 
size [(i?)], the thermal fluctuations Xr — (-^^) ~ (^)^' 
and the disorder fluctuations Xd — [{R)^] ^ [{R)]'^- 

In Fig. JTJJa), we plot the average domain size [{R)] 
as function of the pinning concentration p. As might 
be expected, the domain size decreases with increasing 
concentration. Following theoretical predictions for the 
random-flcld Ising model [S7H5^ . we anticipate that 



(28) 



In [{R)] ex p- 



with exponent p > 0, whereby we assume that [{R}] is 
the analogue of the "break-up" length. For bimodal and 
Gaussian random-fields p = 2, but since it is a priori 
unclear how pinning sites compare to such fields, we leave 
p as a free parameter. The curve in Fig. I7|a) shows the 
corresponding fit of Eq. (28) to our data, where p ~ 0.25 



was used, and the agreement is quite reasonable. 

Next, we consider the magnitude of the thermal xt 
and disorder fiuctuations xd pig-j^^)]- ^^ P increases, 
the fluctuations become smaller, but we always flnd that 
Xd > XT- The disorder fluctuations thus dominate, as 
we had already announced previously. In Fig. mc), we 
plot the distribution P{{R)) of the thermally averaged 
domain size (R) between samples for p = 0.05 (the width 
of this distribution thus reflects the disorder fluctuation 
Xd)- Experiments have indicated that the tail of the dis- 
tribution is exponential PJ. The dashed line in Fig.l7|[c) 
shows a fit to the tail of the simulated distribution using 
an exponential function of the form P{x) ex e"**^. The 
fit captures the data rather well, where 5 ~ 2.1 was used. 

Last but not least, we show in Fig. JT^d) the variation 
of the average domain size [{R)] with the lateral correla- 
tion length ^11 . As ^u increases, there is a mild decrease 
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FIG. 8: (Color online) Simulation results for the pinned 
membrane using pinning height hp = hp — —a/2. In this 
case, there is no random- field effect, (a) The "running av- 
erage" of the spatial fluctuation xs- Clearly visible is that 
Xs — >■ 0. (b) The "running average" of the chemical potential 
fi, which approaches e^ (horizontal line). All data refer to 
p = 0.05, a — 2.0, L — 40. In each plot, results from five 
independent simulation runs are shown. 
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FIG. 7: (Color online) (a) Average domain size [{R}] as 
function of pinning density p (symbols are simulation data; 
the curve is a fit to Eq. ( 28 1). (b) The disorder fiuctuation xd 
and the thermal fiuctuation xt in the domain size as function 
of p. (c) Distribution of the thermally averaged domain size 
between samples for p — 0.05. Note the logarithmic vertical 
scale. The dashed line shows the result of an exponential fit 
(see details in text), (d) Variation of the typical domain size 
[{R}] with the lateral correlation length f y (dots). Also shown 
for each measurement are the magnitudes of the disorder xd 
and thermal fiuctuations xt- (Note: All averages in this fig- 
ure were computed using Eq. (27 1 with T ~ 500, K — 20 
(a,b,d), and K — 500 (c). The data in (a-c) were obtained for 
K = 1, while in (d) k was varied.) 



of the domain size. However, on the scale of the (domi- 
nating) disorder fluctuations, the effect may be difficult 
to observe in experiments. 



Xs ill the thermally averaged bond occupation variables 
(Eq. (24|). The figure strikingly shows xs -^ 0, which 



confirms that the local chemical potential excess has in- 
deed vanished. As for the unpinned membrane, xs now 
reflects the statistical error of our MC simulation, and 
therefore decays ex r~^/^, r being the number of MC 
moves. This result is to be contrasted with Fig. Isja), 
where xs for the pinned membrane converged to a finite 
value. Hence, even though the pinning concentration is 
the same in both cases, choosing hp = h^ completely 
destroys the random-field effect. In Fig. l8[b), we show 
the "running average" of the chemical potential obtained 
via Widom insertion [40l [41] . Interestingly, fj, converges 
exactly onto e™ = — a^/2 of the unpinned membrane. 
Apparently, for neutral pinning sites, the symmetry line 
of the unpinned membrane is restored. 



4. The case of "neutral" pinning sites 

As can be seen in Fig. [5jb,c), for hp = the pin- 
ning sites "repel" closed bonds. When the pinning height 
is sufficiently negative, the reverse situation is obtained 
and pinning sites attract closed bonds (we have checked 
this case for hp = —a). Therefore, it is plausible that 
for some intermediate height hp the pinning sites be- 
come neutral. This special height is given by the sym- 
metry height of the unpinned membrane: hp = —a/2 
(Section IVB2|. Consequently, when hp = hp, we no 



longer expect the local chemical potential excess to be 
a quenched random variable, but rather Afii = for all 
sites. 

To test this assertion, we performed a canonical sim- 
ulation {(/) = 1/2) using pinning concentration p = 0.05, 
pinning height hp = hp, and a = 2.0. In Fig. Isl^a), 
we plot the "running average" of the spatial fluctuation 



At hp — hp, the Imry-Ma argument thus no longer 
applies. We therefore expect a critical point again, at 
a = a' , and macroscopic adhesion domains for a > a'. 
The value of a' should increase with p, and limp_yo a' = 
a* of the unpinned membrane. The quenched disorder 
induced by neutral pinning sites is of the dilution type. 
The diluted Ising model is similar to the standard Ising 
model, but with a (small) fraction of randomly selected 
sites removed from the lattice [42j . Note that this type of 
quenched disorder does not break spin reversal symmetry, 
consistent with our observation that /i -> e^ [Fig.lsjb)]. 

In realistic situations, we do not expect the pinning 
sites to be neutral, nor that the pinning height is the same 
for all pinning sites (as assumed in the present study). 
In all these cases, it is the random-field scenario that ap- 
plies. Nevertheless, for our fundamental understanding 
we felt inclined to discuss the neutral case also. 
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We have studied a minimal model [HES] describing the 
specific adhesion of a cell (or biomimetic vesicle) to an- 
other cell, the extracellular matrix, or a substrate. The 
model incorporates the arguably most important mech- 
anism: the coupling of thermal fluctuations of the cell 
membrane to the state of the receptor-ligand pairs re- 
sponsible for the cell attachment. The latter are effec- 
tively described as either open or closed bonds. Inte- 
grating out the membrane fluctuations this model can 
be mapped exactly onto the two-dimensional lattice gas 
(or Ising model). This mapping confirms previous nu- 
merical evidence of a phase separation [S] between an 
unbound and a bound state. For fixed concentration of 
mobile bonds this scenario implies the formation of a sin- 
gle macroscopic domain of bonds. Moreover, we have ex- 
actly determined the critical point, and we have demon- 
strated numerically that the full model indeed belongs to 
the Ising universality class as anticipated from the theo- 
retical mapping. 

A striking observation in experiments is the absence of 
a single large domain. Rather, adhesion domains of finite 
size form. This observation is often explained as either 
caused by active processes, or the "corralling" of adhesion 
proteins in compartments f21|. This putative hindrance 
of protein diffusion is a purely dynamical approach, which 
does not alter equilibrium properties. Quite in contrast, 
here we have demonstrated that already in equilibrium fi- 
nite sized adhesion domains are implied in a wide variety 
of contexts if one takes into account "pinning sites" that 
locally suppress membrane height fluctuations. These 
pinning sites model perturbations that necessarily occur 
in vivo due to the crowded, highly non-ideal composition 
of cell membranes and the anchoring of the cytoskeleton 
to the ECM or other cells. 

We have presented both analytical calculations and nu- 
merical evidence that the presence of pinning sites corre- 
sponds to quenched disorder which, in the language of the 
Ising model, induces a random field that prevents macro- 
scopic domain formation in two dimensions ^ 2^25] . This 
is in sharp contrast to the type of quenched disorder 
where the positions of ligand-receptor bonds are ran- 
dom. The latter constitutes a marginal perturbation [11] 
and does not fundamentally alter the scenario of Fig. |4j 
i.e., macroscopic domain formation above some critical a. 
In contrast, random- field disorder is a relevant perturba- 
tion. It requires that the disorder couples linearly to 
the order parameter. We have demonstrated here that 
this condition is typically fulfilled for membrane adhe- 
sion, where membrane height pinning (disorder) couples 
to bond formation (order parameter (p). Hence, the fun- 
damentals of statistical physics can be used to explain 
adhesion domains of finite size in equilibrium. 
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Appendix A: Height correlations 

Using periodic boundaries the separation profile /i(r) 
is expanded into Fourier modes through 



Mr) = VV"'''", 



q 



hq^jf d^r h{r)e''^' 



with hq = h' + ih'' Since the separation field is real 

h{r) = ho + 2^[h'^ cos(q ■ r) + h'^ sin(q • r)], 
qeQ 

where the set Q contains the independent wave vectors 
excluding q — 0. 

For a free membrane in the absence of both bonds and 
pinning sites the Hamiltonian Eq. (fl]) in Fourier space 
reads 

qGS 

with 7o = PA-y and 7^ = 27o[l -I- (9^11)"*]. Hence, the 
mean height is (/i(r))o = and the height correlations 
read 

n n-. /,/N,/A\ 1 , v^ cosq- (r- r') 
m r - r' = /i r /i r' = - + 4 ^J ^"^ -■ 



70 



qeQ 



The height correlations between two points on the mem- 
brane decay on the length scale ^y . For large separations 
they reach m{r — > 00) = (/io)o ~ I/70 due to the zero 
mode fluctuations. For large A we can neglect this con- 
tribution. Replacing the sum over discrete wave vectors 
by an integral we then obtain 



m{r) 
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qJoiir) 



2n(inJ^ M + (<7e||)4' 



where Jq is the zero-order Bessel function of the first 
kind. Performing the final integration leads to Eq. (pi) 
for the height correlations of a free membrane. 



Appendix B: Effective free energy 

The full partition sum of the system reads 

• • • atiN e"''^" 



Z = ^ f[dh{r)] f dhi---dh 



N 



M 



l[S{h{r,)-K) II £S{h{rk)-hp), (Bl) 

i = l k=N+l 
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where the Hamiltonian "Hat is given in Eq. (Is]). The 
M = N + n total constraints due to the N bonds and 
the presence of n pinning sites are represented through 
(5-functions, where ^ is a microscopic length. We normal- 
ize the functional measure such that Z ~ I for a free 
membrane. 

As usual, we represent the J-functions as integrals 



S{hir,) - hi) = 



1 

2^ 



dA,, 



,iA.[/i(r.)-M 



where we complement hi — hp for i > N . This allows 



us to write the partition sum Eq. (Bl I involving a single 
exponential 



Z = Y, [Ah{v) 
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with coefficients 
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Co = i^Ai - /3a^bi, 

1=1 1=1 

N 

Cq = 2iyj Aicosq • r^ — 2/3a V] ^i cos q • r^, 

i=l i=l 

A/ AT 

Cq = 2i 2^ Ai sin q • r^ — 2/3a ^^ ^i sin q • r^. 



M 



Note that the first sum runs over all sites whereas the sec- 
ond sum only includes the sites where bonds are present. 
Performing the Gaussian integrations over the height 
modes {/iq} and the auxiliary variables {Xi} we obtain 

fn r 

^ ^ (2^)M/2 ^dct m-1 Y^jdh,--- Ah^ 

\ I ^' ^ 1 

X exp<f "2 X^^"^ ^)ij^3^3 - /3^bi[ahi - £b] > • 

(B2) 
The prefactor is independent of the state of the bonds. 



qeQ 






N 



While in principle it depends on the quenched disorder 
of the pinning sites it does not contribute to thermal 
averages. Setting hp — the terms with indices i > N 
corresponding to the pinning sites drop out of the sum 
and we obtain Eq. Q. 



Appendix C: Mean height 

The derivation presented in the previous section is 
easily modified to calculate a generating function, from 
which arbitrary moments can be obtained. For example, 
replacing cq — ?> cq -I- A we obtain the generating function 
Z(A) from which the mean height follows as 



(ho) 



ainZ(A) 



dX 



(CI) 



A=0 



Repeating the calculation for n = we obtain 

{ho)=-^{^). 



(C2) 
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